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SUMMARY 


A method is presented which relates the aerodynamic damping derivatives 
as a function of instantaneous angular displacement with the aerodynamic deriv- 
atives as a function of oscillation amplitude. Each derivative is expressed 
in a power series, and By applying the Rryloff-Bogoliuboff equivalent linear- 
ization technique, the solution is obtained by equating like coefficients. It 
is assumed that the damping per cycle is small, the damping coefficient is 
symmetrical, and the pitching moment is linear. Comparisons with numerical 
integration results show no appreciable error. 

The effect of a nonlinear pitching moment is considered for a represent- 
ative case and is found to be small except when approaching unstable trim. 


INTRODUCTION 


The aerodynamic pitch damping coefficient can be obtained from wind-tunnel 

tests by two general methods. In one method, an example of which is small- 

amplitude forced-oscillation tests (ref. l), the aerodynamic pitch damping 

coefficient C + C is measured as a function of instantaneous angular 
m m. B 

q a 

displacement 9 by forced oscillations of small amplitudes at discreet angles 

of attack. In the other method, examples of which are ballistic, free-flight, 

or free- oscillation wind-tunnel tests (ref. 2), the average aerodynamic pitch 

damping coefficient C— + C— is measured as a function of oscillation peak 

in m. 

q a 

amplitude 9 q by determining the change in peak amplitude per cycle. 

To be used in current digital trajectory computer programs the aerodynamic 
pitch damping coefficient must be expressed as a function of instantaneous 
angular displacement rather than oscillation peak amplitude. Since ballistic. 


free -flight, and free- oscillation wind-tunnel experiments determine the damping 
as a function of oscillation peak amplitude, a method of expressing these re- 
sults as a function of instantaneous angle of attack is needed. It is the 
purpose of this paper to develop a method whereby damping as a function of 
oscillation amplitude can be expressed as a function of instantaneous angle of 
attack. The method has the added capability of determining the damping as a 
function of oscillation amplitude when damping as a function of instantaneous 
angle of attack is known. This latter property allows an oscillation-amplitude 
time history to be constructed without resorting to numerical integration of 
the equation of motion. 


SYMBOLS 


A, Ag, A^, Ag, Aq 

B, B 2 

c 2 , C 4 , Cg, C 8 
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f(e) 
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q 


arbitrary constants in the series expansion of f ( 9 ) 
arbitrary constants in the cubic series for C 

m 

arbitrary constants in the series expansion of f^0^ 

pitching-moment coefficient, — -^ c . ^^ n S -. ™ amen '^ 

dC 

m 

ctec 

/ d C \ / dC \ 

damping in pitch coefficient, — 2 + _J£ 

\^v/ \^v/ 

average damping in pitch coefficient over one full 
oscillation 

reference length 

viscous damping as a function of 9 

viscous damping as a function of 9^ 

mass moment of inertia 
slope of assumed linear pitching moment 
power per cycle as a function of f(9) 

power per cycle as a function of 
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angular pitching velocity 
dynamic pressure 


s 

T 

t 
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lc 

0 * 

lc 


ut 
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reference area 

period 

time 

free- stream velocity 
angle of attack 

rate of change of angle of attack 
instantaneous angular displacement 
oscillation peak amplitude 

limit- cycle amplitude, linear pitching moment 
limit-cycle amplitude, cubic pitching moment 
unstable trim amplitude, cubic pitching moment 
angular velocity 
angular acceleration 
arbitrary angular displacement 

2tc 

angular frequency, 


ANALYSIS 


Derivation of Equations 

The second- order differential equation of motion for a single-degree- of - 
freedom oscillating system with viscous damping proportional to displacement is 

10 + f(e)e + ke = o (i) 

where f(9) is the instantaneous viscous damping and k is the slope of an 
assumed linear pitching moment. The function f(9) corresponds to the damping 
coefficient as found in a forced oscillation experiment and, for a vehicle 
symmetrical in the pitch plane, can be expressed as 

f(9) = A + A ^ 9^ + A^9^ + Ag0^ + Ag9^ + . . . (2) 
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If, instead of small-amplitude oscillations about a fixed mean angle of 
attack, the vehicle had oscillated harmonically about the trim angle with a 
large amplitude, the power per cycle would be 


p = f f(e)e 2 dt - / l a + A 2 e 2 + A 4 e 4 + A g 

Jo Jo ' 

For a free oscillation, the equation of motion is 

10 + f(e 0 )e + k 0 = 0 


0 6 + Ag0 8 + . . .^ 9 2 dt 




( 3 ) 


(4) 


where is " fche viscous damping as a function of oscillation amplitude. 

The function ^(®o) corresponds to the average damping coefficient as deter- 
mined in a free-oscillation experiment, and for a vehicle symmetrical in the 
pitch plane, can be expressed as 


f(e 0 ) . o + * o u eJ + Cg9g + o 8 e® + 

The pcwer per cycle can be written as 


( 5 ) 


? 0 ■ f ( 6 o) - ( 0 + Vo 


+ Vo + °6 S 0 + °8 e 0 + 


*'o 


2 

e dt 


(6) 


Kryloff and Bogoliuboff (ref. 3) have shown that equation ( b ) is the equivalent 
linear differential equation of the nonlinear equation of motion, equation (l). 
As in reference b, where the power per cycle is assumed to be the same over a 
given oscillation, then 


Solving for f ( e 0 ) 


f f(e)e 2 dt = f (e 0 ) f i 

J o J o 


p 

e^dt 


(7) 


/ f(e) 0 ‘ 
f ( e o)= ° 


'dt 


f T. 

J e 

J 0 


e 2 dt 


( 8 ) 


b 


or its equivalent 


n 


c + c 2 e o + Vo + c 6 0 o + c 8 e 0 + • • • = ~ 


For almost harmonic motion 0 can he replaced hy 


A + A^9 2 + A^0 + Ag9 + Aq 0 + . . . ^ 9 dt 


a 

Jc\ 


e 2 dt 


(9) 


9 = 0 Q cos at 


(10) 


If the damping moment is small, then 9 will he a slowly varying function 
of time during the period T and as a first approximation can he considered 
constant over a cycle. The relation of equation (10) allows a straightforward 
integration of equation (9)* Integrating and carrying out the indicated divi- 
sion 


c + C 2 e o + <Vo + °6 e o + °8 e o 

Equating like coefficients 


A + sVo + K e o + 


l 6 o 


7 . 


+ lfeV 


C 2=K 

C ll " 5% 

c 6 

°8 ' 325*8 

Terms heyond the eighth order were not considered. 


(ID 


( 12 ) 


If, from a forced-oscillation test, f(9) is known in a power series such 
as equation (2), then can ^ e ^ erm ^- ne( ^ directly as 

f ( 9 o) = A + ¥ A 2 e 0 + SVo + 0 + 128 A 8 9 0 (15) 

On the other hand, if f(®o) known from a free-oscillation test and can he 
expressed in a power series such as equation ( 5 ), f( 0) can he found directly 
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as 


f(e) = c + 4 c 2 e 2 + 8 c 4 e 4 + yc 6 e 6 + ^ 8 e 8 


f(e) = 


(c + c N 

m m. J 

V <1 a/ 


' ( 6 o) = ~ 2 V f C m + C m. 


' V <1 a,/ J 

the desired relationship between C + C and C_ + C_ is found. 

mm. mm. 

q a q a 


Methods of Application 

Determination of C_ + C_ from measured values of C + C .- If 

in m. m m. 

2, a q a 

values of C + C as a function of angular displacement 9 have heen 
m m. 

q a 

measured in small-amplitude forced-oscillation tests, the damping coefficient 
can he expressed in a power series 

c + c - gv (a + AJ 2 + + kA + Ac.e 8 ') 

"<1 m a ^SD 2 V 2 ** 6 8 1 

Then, from equation (13) the damping coefficient as a function of amplitude 
of oscillation is 


C m + C m. = ( A + iVo + BVo + ^ A 6 9 0 + lfeV 

q a q^SD \ 


With this calculated C_ + C_ the oscillation-amplitude time history can he 

in m. 

q a 

determined hy using the logarithmic decrement relation 

^ o) - (\ + % ) (16' 
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< 0) 

= 0 , 

0,2 

0,1 


0,1 

If the measured 

values of 

C 

m 

+ C 

m. 


q a 

amplitudes and positive damping at larger amplitudes, a limit cycle is indicated. 
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The amplitude of this limit cycle is the amplitude at which the oscillation- 
amplitude time history, determined "by using equation (l 6 ), levels out. The 
amplitude of the limit cycle can also he found by solving for the roots of 
the power series representation of C_ + C_ . The two real roots, opposite 

q d 


in sign but equal in magnitude, near trim will be the amplitude of the limit 
cycle . 


Determination of 


from measured values of C_ + C_ . 

m m. 

q a 


If an 


oscillation-amplitude time history such as that shown in sketch (a) is obtained 
from a free- oscillation test, the damping as a function of oscillation ampli- 
tude can be calculated with the following equation: 


’o,2h - 




Sketch (a) 

Then, ^(®q) can ”' :>e Pitted against 9 ^ and expressed in the power series of 

equation (5) • Damping as a function of angular displacement can be calculated 
from equation (l^). Prom equation (15), the damping coefficient as a function 
of angular displacement is 


Numerical Comparisons 

The following two theoretical methods have been presented to relate the 
pitch damping derivatives as a function of angular displacement with the pitch 
damping derivatives as a function of oscillation amplitude: 


(1) 

Determination 

of 

C_ + 
m 

q. 

C- 

m. 

a 

from measured 

values 

of 

C + 
m 

q 

c 

m. 

a 

(2) 

Determination 

of 

c + 
m q 

C 

ni- 

ce 

from measured 

values 

of 

c_ + 

m 

q. 

CL. 

m. 

a 


In order to ascertain the accuracy of these methods, arbitrary damping deriv- 
ative coefficients C + C were programed in the single-degree-of- freedom 

mm. 
q a 

equation of motion 

10 - — ( C „ ♦ V - ° = ° (l8) 

V <1 a/ a 

and numerically integrated. The oscillation-amplitude time history obtained 

was used to calculate C_ + C_ . The results of these calculations were 

m m. 

q a 

then compared with the theoretical solution to verify method 1. To ascertain 

the validity of method 2, C + C was obtained from the calculated values 

m. m. 

q a 

of C_ + C_ and then were compared with the input values, 
q m d 


The following two cases present 
numerical comparisons made to verify 
the two methods. 

Case A . - The first case is that of 

a quadratic with positive damping. The 

value of C + C was chosen as 
m m. 
q a 

-6. 55 X 10-5 Th e oscillation- 

amplitude time history in figure 1 was 

obtained by numerical integration of 

equation ( l8) - and from the time history, 

C_ + C_ was calculated by using equa- 
m m. 
q a 

tion (17). The theoretical damping co- 
efficient was then determined by using 
method 1. The calculated and theoretical 
values are compared in figure 2. 

Method 2 was used to obtain 

C . C from the calculated values of 
in + m. 

1 a 

C_ + C_ of figure 2, and the results 
mm* 
q a 

are compared with the actual value 
-6.55 x 10” 5 9^ in figure 3- 
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Figure 1. - Oscillation-amplitude time history 

for case A, C m + C m« = “** 55 X 10 3 ® 2 ‘ 
q « 
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Case B . - The second case is that of a quadratic vith both positive and neg- 

mm 'z o 

ative damping. The value of C + C was chosen as 2.0 - 6.55 x 10 J 9 

mm. 
q. a 

Equation (l8) was numerically integrated to obtain the oscillation-amplitude time 

history of figure 4 which shows a limit cycle at 55° • From figure 4. C_ + C_ 

m m. 

q a 

was calculated by using equation (l6) , and the results are compared in figure 5 
with the theoretical coefficients obtained from method 1. Note that both the 
theoretical and computed values show zero damping at 9^ = 35% as would be 
expected. 

Method 2 was used to obtain C + C from the calculated values of 

m m. 
q a 

C_ + C_ , and the results are compared with the actual coefficient 
in Iff. 

q a 

2. 0 - 6.55 * 10 ^ 6 in figure 6. Both solutions have damping coefficients of 
0 at 9 = 17. 5 °, which is exactly one-half the amplitude of the limit cycle. 
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Figure 4. • Oscillation-amplitude time history lor case B, C m + C m = 2 -6.55 * 1O - ^0^ 

1 « 
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One of the important applications of method 1 is the construction of an 

oscillation-amplitude time history when C + C is known hy numerical inte- 

hi ™ ^ 


m. 

a 


gration of equation (l6) rather than numerically integrating the entire equation 


hy method 1. Then, 


of motion. First, C + C is converted to C_ + C_ 

m m* m m. 

q a q a 

hy using equation (l6), the oscillation time history can he constructed with 

small time steps so that both 0 and C_ + C_ are considered to he constant. 

u in m. 

q a 

This application is illustrated for case B where: 


C + C = 2.0 - 6.55 x 10-5 e 2 
in m. 

q a 


By using method 1 


C S + c m. - 2 '° -(y) 6-55X10-5 e 


q. a 

This is shown in figure 5« Substituting this value in equation (l6) 

-.2 


6 0,2 6 0,1 e 


§ izr- [ 2 - (s) 6 -55 x 


The time step was altered so that 9 0 changed approximately 1° each step, and 

the oscillation-amplitude time history obtained is shown in figure 7. The 
small differences between the actual and constructed time histories are attrib- 
uted to the assumption that C_ + C_ is constant over each time step. A 

m n 

q ex 

limit cycle of 55° is also shown in this figure. The amplitude of the limit 
cycle can also be found by solving for the roots of C_ + C_ 


\ + \ - 1ST- [ 2 - (I) 6 -55 x 10' 3 = 0 

2 2 
0 5(6.55 x 10" 5 ) 

e 0 - * 55 “ 
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Effect of a nonlinear pitching moment .- In the development of the theory 
in this paper, only linear pitching moments were considered. However, except 
for small oscillations, the pitching moment for many vehicles is nonlinear with 
amplitude. To determine the usefulness of this theory as applied to the non- 
linear pitching moment, the amplitude of a predicted limit cycle was compared 
with an actual limit cycle, since this was believed to be the most severe lim- 
itation. The nonlinear pitching-moment coefficient was simulated in the single 
degree- of- freedom computer program by 

C = B0 - B e 5 
m 2 

which is shown in sketch (b) where unstable trim is denoted as 


e ut " ± 


15 





as shown in sketch (d) where ±2 ~\J A/A^ denotes the amplitude of the theo- 
retical limit cycle 9^ . Various values of A were put into the computer 
program to vary the amplitude of the actual limit cycle 0“ from 0 to the 

ft 

point at which tumbling occurred. The values 0.^ and 0^ c were nondimen- 
sionalized by dividing by 0 


P a 



From figure 8, it can be assumed that the nonlinear pitching moment would 
cause no appreciable error in the solution for the amplitude of the limit cycle 
up to values of 0.7 \ c /® u -f Fran 0.7 to 0*9 ®i c /® u t> 'the error 



Figure 8. - Comparison of actual limit cycle with theoretical limit cycle. 
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Unstable 



increases; and for values above 0. 9 the linear approximation is not 

usable since it predicts a limit cycle when, actually, tumbling occurs. 


CONCLUDING REMARKS 


The aerodynamic pitch damping coefficient as a function of instantaneous 
angular displacement is related with the aerodynamic pitch damping coefficient 
as a function of oscillation amplitude by a method in which the Kryloff- 
Bosoliuboff equivalent linearization technique is used. Power series expres- 
sions for the functions are developed and the like coefficients of the series 
are equated to provide the relationship. It is assumed that the damping per 
cycle is small, the damping coefficient is symmetrical, and the pitching moment 
is linear. Comparisons with numerical integration results show no appreciable 
error. The effect of a nonlinear pitching moment is considered for a repre- 
sentative case and is found to be small, except when approaching unstable trim. 
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